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Abstract 

Rotation of tokamak-plasmas, not at the mechanical equilibrium, is investigated utilizing the 
Prigogine thermodynamic theorem. This theorem establishes that, for systems confined in rect- 
angular boxes, the global motion of the system with barycentric velocity does not contribute to 
dissipation. This result, suitably applied to toroidally confined plasmas, suggests that the global 
barycentric rotations of the plasma, in the toroidal and poloidal directions, are pure reversible 
processes. In case of negligible viscosity and by supposing the validity of the balance equation 
for the internal forces, we show that the plasma, even not in the mechanical equilibrium, may 
freely rotate in the toroidal direction with an angular frequency, which may be higher than the 
neoclassical estimation. In addition, its toroidal rotation may cause the plasma to rotate globally 
in the poloidal direction at a speed faster than the expression found by the neoclassical theory. 
The eventual configuration is attained when the toroidal and poloidal angular frequencies reaches 
the values that minimize dissipation. 
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I. INTRODUCTION 



Toroidal and poloidal plasma rotations are of great importance in fusion science as it plays 
a very important role for controlling transport in magnetically confined tokamak-plasmas 
PP. Apart from improving the confinement, intrinsic plasma rotation is also important be- 
cause it may affect the threshold to an improved confined state (L-H transition) [2] and, if 
combined with large rotational shear, it can produce the formation of an internal transport 
barriers (ITB) by suppression of turbulence [3]. In addition, rotation is of significant interest 
for ITER performance. Recent experiments in the DIII-D tokamak have indeed shown that 
plasma rotation may suppress the resistive wall mode (RWM) in tokamaks |4j. The rota- 
tional velocity observed in tokamaks may be very high, especially in the toroidal direction, 
which may reach values of the speed of sound. Rotation in the poloidal direction is however 
much slower. In tokamaks, torque yielding rotation can have several origins, among which, 
heating systems such as Neutral Beam Injection (NBI) (see, for example, Ref. [5]) and radio- 
frequency heating (see, for example, Ref. [6]). However, even in absence of external torques, 
axisymmetric toroidal plasmas exhibit an intrinsic rotation both in the toroidal and the 
poloidal directions [7]. Several mechanisms have been proposed for explaining the origin of 
the intrinsic plasma rotation. Without claim of completeness, we cite three of them. We 
first mention the theory where the turbulence driven residual stress, acting as an intrinsic 
torque, is supposed to spin-up toroidal rotation [8] (see also Ref. [9] for an overview). In 
this approach, the toroidal moment transport is driven by the parallel and the perpendicular 
Reynolds Stress, -< u r U0 y (with u r and u^ denoting the barycentric velocity in the radial 
and the toroidal directions, respectively, and -<••■>- the mean value operation) [TU] . Ac- 
cording to this theory, the non-diffusive residual stress is responsible for the intrinsic plasma 
rotation, through a ky -symmetry breaking mechanism able to convert radial inhomogeneity 
into parallel spectral asymmetry [11]. Another approach can be found in Ref. [12]. This 
theory is based on the transfer of angular momentum to the wall and towards the centre of 
the plasma column by modes associated with the gradients of the plasma temperature and 
density. The inward transport of angular momentum from the edge of the plasma column to 
its centre is accounted for by the modes that are associated with the gradients of the plasma 
temperature and density in the presence of a velocity flow along the magnetic field. Finally, 
we cite the mechanism reported in Ref [13]. In this theory, the resultant plasma toroidal 
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rotation equation includes the effects of collision-induced perpendicular viscosities, anoma- 
lous transport due to micro-turbulence in the plasma and momentum sources. According to 
this picture, non-axisymmetric fields produce a toroidal torque throughout the plasma that 
relaxes the toroidal flow to an intrinsic ion-temperature-gradient diamagnetic-type flow in 
the direction counter to the plasma current. 

In this work, we do not address the validity of the (several) mechanisms proposed in literature 
for explaining the onset of the plasma spin-up. This paper treats (and solves) the following 
crucial questions on Tokamak-plasma rotation in absence of external torque. 

Even not in the mechanical equilibrium, 

- Why does the plasma may rotate without dissipations ? 

- Why does the frequency in the toroidal direction may largely be different from the neoclas- 
sical estimation ? In this case, what is the role of the initial conditions ? 

- Why, on the contrary, the value of the frequency in the poloidal direction tends to be damped 
towards the neoclassical value ? 

To the best of our knowledge, we cannot find in literature a rigorous answer to this series of 
questions. It should be noted that our results derive directly from rigorous thermodynamic 
theorems and they are not based upon a theory. In particular, our scope is to illustrate the 
Prigogine macroscopic thermodynamic theorem (PMT) which, in our opinion, may signifi- 
cantly contribute to the understanding of the features of the intrinsic rotation of toroidally 
confined plasmas. The PMT asserts that, for systems confined in rectangular boxes, the 
global motion of the system with barycentric velocity does not contribute to dissipation 
[22], [IS]- We shall adopt the single fluid picture and focus on the ion fluid, whose velocity 
can be considered equal to that of the single fluid, because the electron mass is negligible 
compared to the ion mass. By combining the results reported in Refs [IE] and p2] with 
the PMT, suitably readapted to tokamak-plasmas, we show that the barycentric rotation of 
dissipative tokamak-plasmas is a pure reversible process. In addition, if viscosity is negligi- 
ble, we demonstrate that the global motion of the magnetically confined collisional plasmas, 
satisfying the balance equation for the internal forces, is determined up to an arbitrary flow 
velocity. In this respect, to avoid misunderstanding, we would like to clarify the following. 
It is known that an immediate application of the ideal magnetohydrodynamics to plasmas, 
at the mechanical equilibrium, leads to the result that the barycentric motion of the system 
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is determined up to an arbitrary velocity. However, here we prove a stronger result. In 
absence of viscosity, the global motion of collisional plasmas is determined up to an arbi- 
trary velocity, even if the system is not at the mechanical equilibrium. This theorem is a 
kind of thermodynamic equivalency version of the mechanical inertial principle, according 
to the scheme: Inertial Principle f± Reversible Process. A direct consequence of this theo- 
rem, applied to tokamak-plasmas, leads to prove that the plasma may freely rotate with a 
toroidal angular frequency, which may be higher than the neoclassical prediction, with an 
arbitrary additional factor equal to Aus^ ~ v^/R (with v^ and R denoting and arbitrary 
toroidal velocity and the major radius of the tokamak, respectively). In its turn, the toroidal 
rotation may cause the plasma to rotate globally in the poloidal direction at a speed faster 
than the expression found by the neoclassical theory, with an additional contribution equal 
to Aug = V(/,Bo/ \rBij,) (where Bg and B^ stand respectively for the poloidal and the toroidal 
magnetic field configuration, and r is the radial coordinate). Let us mention that the ther- 
modynamic approach explains the Solomon experiment, without invoking the existence of an 
"anomalous torque" [TS]. Indeed, Solomon et at, by studying the momentum confinement 
for the DIII-D tokamak as a function of applied neutral beam, found that, under balanced 
neutral beam injection {i.e. under zero total torque), the plasma still maintains a significant 
rotation in the co-direction. A preliminary work aiming to understand the plasma rotation 
effect through a thermodynamic approach, can be found in the footnote [37] . 
The paper is organized as follows. In Section [IT] we discuss the PMT and its main conse- 
quences. In particular, the Subsection II A| deals with dissipative systems, which are not 



in mechanical equilibrium. Conversely, the Subsection |II B| refers specifically to conserva- 

is 



tive diffusion processes occurring in systems at the mechanical equilibrium. Section III 



devoted to the application of the MPT to the intrinsic rotation effect in tokamak-plasmas. 
The frequency rotation in the poloidal direction is obtained in the Velikhov geometry [19] 
by removing any supplementary ad hoc hypothesis. Concluding remarks can be found in 
Section HVl 
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II. BASIC EQUATIONS FOR A SINGLE FLUID DESCRIPTION AND PRI- 
GOGINE'S THEOREM 



We adopt the single fluid description, which excludes the issue of the momentum transfer 
between the ion and electron fluid. In order to set up the nomenclature, we briefly recall 
the basic definitions below. We assume that the barycentric velocity satisfies the evolution 
equation [20J 

,-v.n + * (i) 

with II denoting the pressure tensor measured by an observer moving with the barycentric 
velocity u j3H]- Also the substantial time-derivative, d t , is constructed with the barycentric 
velocity : 

The single fluid quantities, such as the barycentric velocity, u, the mass density, p, the 
barycentric force, J 7 , and the (density of ) flux diffusion of species a, are defined in the 
usual way 

777-jUj -|- Tfl e Xl e v -\ v -\ v -\ p a , . 

u = ■ ; p = ) p a = } n a m a ; 7 = ) — T a (3) 

rrii + m e p 

a=e,i a=e,i a=e,i 

Ji = p a (u a -u) ; a = (e,i) 

Here, m 8 and m e , and Uj and u e are the ion and the electron mass, and the ion and the 
electron velocities, respectively. q e denotes the absolute value of the electron charge and n a 
is the species number density. For the quasi neutrality condition, we have n e = Zui ~ Zn, 
with Z denoting the charge number. 



A. Dissipative systems not in mechanical equilibrium 



The PMT holds by assuming the validity of the following hypothesis : 

The Gibbs relation, valid for closed and open systems, also applies to the barycentric motion 
of the single fluid, under the condition that the derivatives appearing in this equation are 
modified according to the motion of the fluid i.e., 
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ds du dv ^-^ dN a 

T dt = ti +p Tt-^^ (4) 

a=e,i 

with the total derivatives provided by formula §2^. 

s and n a indicate the total entropy of the system per unit mass and the chemical potential 
(for species a) per unit mass, respectively. T, u, v and N a are the temperature, the energy 
density per unit mass, the specific volume (v = p _1 ) and the mass fractions, respectively. 
Prigogine demonstrated that, under the validity of this hypothesis, the entropy production 
a, related to a continuos media, can be brought into the form [22] 

d fl\ 1 x - du K 1 v - x - r . fdfi c 



a 



j=l,2,3 J j'k=1,2,3 J a=e,i j=l,2,3 J 



"4EE[^-(|) r -^ <«> 

a=e,i j =1,2, 3 J 

Here, indexes j, k distinguish the components of a vector or of the viscous tensor 7T. Moreover 
(xi,x 2 ,x 3 ) stands for (x,y,z) [i.e., the three spacial cartesian coordinates]. R a is the total 
force (with components B? a ) exerted on the particle species in question as a result of collision 
with other species in the plasma and u K denotes the components of the barycentric velocity. 
Hence, we have R Q = J2p where R a/ 3 represents the interaction between species a and 
P. 3 h denotes the entaltpy-heat flow linked to diffusion defined as 

j h = r + (6) 

a=e,i 

with h a denoting the enthalpy per unit mass of species a, 3 q the heat flow, and the flux 
diffusion of species a, J^, is defined in Eqs respectively. The first term on the right hand 
side of the first equation of Eqs ^ expresses the dissipation due to the heat flow linked to 
diffusion, the second term represents the dissipation by viscosity of the barycentric kinetic 
energy, and the third term takes into account the dissipation by diffusion of the kinetic 
energy. Finally, the last contribution, Tad in Eq. ([5]), is the source of the kinetic energy due 
to diffusion. Notice that the barycentric velocity does not appear explicitly in the entropy 
source strength Eq. (j5|. Expression ({5j contains only the derivatives d Xi vP and the diffusion 
vectors [39]. In other words, the global motion of the system with barycentric velocity 
u does not contribute to the entropy production. Hence, the global motion of the system 
constitutes an "intrinsic" reversible processes. This is the essence of the PMT. This theorem 
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reads as follows. "For systems confined in a rectangular box, global motion of the system 
with barycentric velocity u does not contribute to dissipation" [22], |15j . 
In neoclassical theory, plasma is considered to be in Onsager's region (i.e., the closure re- 
lations are linear or, in other words, the thermodynamic flows are linearly connected with 
the thermodynamic forces) [23]. According to the Minimum Entropy Production theorem 
[Hj, in this region, for time-independent boundary conditions, the state of minimum en- 
tropy production turns out to be a stationary state. Because of Eq. ^ combined with the 
minimum entropy production theorem, also the steady-state reached by the system does not 
depend on the barycentric velocity (but only on its derivatives). Hence, with the addition 
of the Minimum Entropy Production Theorem, the PMT yields that in the Onsager region, 
the steady-state, obtained by minimizing the entropy production, remains unchanged by 
performing a global motion of the system with barycentric velocity. 

We immediately observe that in absence of viscosity and, for collisional plasmas where the 
balance equation for the internal forces is satisfied 



the entropy production of the system, Eq. ([5]), is invariant under the barycentric velocity 
transformation: u — > u + v, with v denoting an arbitrary velocity. Hence for collisional 
plasmas, with no viscosity, global motion of the system is determined up to an arbitrary 
velocity v. Notice that it is not required the condition that the system is at the mechanical 
equilibrium. 

B. Conservative diffusion processes occurring in systems at the mechanical equi- 
librium 

In this subsection we shall consider the particular case of conservative diffusion processes oc- 
curring in systems at the mechanical equilibrium. We briefly check that, in these conditions 
and in case of no viscosity, the balance equation for the internal forces, Eq. ([7]), is automat- 
ically satisfied. The source of the kinetic energy linked to diffusion, Tad, consists of three 
contributions. The first term takes into account the work made by the external forces and 
the second term is the contribution due to the gradients of concentration. Finally, the last 
term represents the dissipation of the diffusion kinetic energy due to the collisions between 




(7) 
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particles with different velocities. In many physical examples, the latter term compensates 
exactly the first two contributions i.e., 

Consequently, ad = 0. In literature, when condition ^ is satisfied, the process is referred 
to as diffusion conservative process. Notice that, in case of diffusion conservative process, 
the entropy production given in Eq. ^ takes the form usually reported in literature (see, 
for example, Refs [15] and [M] ) 

J=l,2,3 J jre=l,2,3 J a=e,i j=l,2,3 J 

When the barycentric acceleration a, in Eq. Q vanishes, the system is at the mechanical 
equilibrium. In absence of viscosity, the mechanical equilibrium condition reads 

£ p°rt = |r ( 10 ) 

a=e,i ^ 

However, from the Gibbs-Duhem relation 

p a dfi a = —psdT + dp (11) 

a=e,i 

we immediately obtain 

£4|;)t = I^ (12) 

a=e,i J J 

Hence, the mechanical equilibrium condition may be rewritten as 



If the diffusion process is conservative, the mechanical equilibrium condition simply reads 



where Eq. (|8j) has been taken into account. Hence, under the above strict restrictions, we 
re-obtain, as a very particular case, the result of the ideal magnetohydrodynamics applied 
to plasmas. At the mechanical equilibrium, global motion of a viscousless fluid, is deter- 
mined up to an arbitrary velocity v. Notice that in this particular case v should satisfy 
the supplementary condition d{v = and its value is determined by specifying the initial 
conditions. 
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III. PLASMA ROTATION 



In this session we shall apply the results discussed in Section [IT] to magnetically confined 
plasmas. First of all, we should check that the PMT, which has been proved for systems 
confined in a rectangular box, applies also to plasmas confined in toroidal geometry. 
• Entropy Production of L-Mode Tokamak-Plasmas 

Let us consider, for simplicity, fully ionized tokamak-plasmas, defined as a collection of 
magnetically confined electrons and positively charged ions, with the magnetic configuration 
written in the form [16] 

B = — x V0 + F(ip)V(f> (15) 

Here ip indicates the poloidal magnetic flux, F is the characteristic of axisymmetric toroidal 
field depending on the surface function, ip, and <p is the poloidal angle, respectively. A 
fundamental remark is the fact that, even in a long mean free path regime, the (density 
of) distribution functions / a (v,x, t) of a magnetically confined plasmas are dominated by a 
local equilibrium reference state and, in the phase space {v,x}, it depend on velocity only 
on through the diffusion vectors (see, for example, Ref. [23] ) 

f (v,x,t) = /°°(JJ,x,t)[l + x"(J«,x,f)] with 

r (js. x, t) = m*. t) ( 2 ^ y f exp - [ maB , (Xi 1 t)ra(Xit) 1 Ji r] (is) 

The deviations from this state, x° > can De expanded in a series of irreducible tensorial 

Hermite polynomials, H[™). iq , getting 

00 00 
X a (3 d a , x, t) = 9° (2k) (x, t)H^\3 d a ) + C (2K+1) (x, t)H^ +1 \3 d a ) (17) 

K=0 K=0 

OO 



+ E<2 K) ( x ' t )<?( J ") + 



with g Q ^ ^ = g a( ^ 2 ' ) = q"^ = and summation over repeated tensor indices, l±,L2, 



is 



understood. The dimensionless entropy production of species a, a a , is derived under the 
sole assumption that the state of the quiescent plasma is not too far from the reference 
local Maxwellian. In the local dynamical triad [25], and up to the second order of the drift 
parameter e = Pl/L <^ 1 [with p L and L = V -1 denoting respectively the Larmor radius 
and the length scale of variation of macroscopic plasma parameter (such as temperature, 
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density, etc.)], a a can be brought into the form (see Refs [16] and [17] ) 

^ = - + cv 3) + ^; (3) ) + w - + w 50 + < 3) ) 
^ = ?£^ (3) + ^; (3) ) + 4 3) y (3) + ^; (3) ) + w° + 

Here gf rc [with the index t taking the three possible directions: radial (p), tangent (||) and 
per-tangent (A)] denote the Hermitian moments of the distribution functions and gt , g?^ 
are the dimensionless source terms. Index k takes the values k = (1, 3, 5). The source terms 
are related to the spatial inhomogeneity, and are expressed as gradients of the fields (i.e., 
Vn Q , VT a , Vu„ etc.). Hence, they depend on space and time (x, t) only. We note that, 
the terms involving the radial thermodynamic forces are identical to those of the classical, 
thermodynamical form given by Eq. ([5]). Hence, in the drift approximation, the entropy 
production is obtained as the sum of two contributions: a classical piece, and a neoclassical 
piece. 

o° = o5 + o% ; a={e,i) (19) 



In the (linear) Onsager region, Eqs ( 18 ) reduce to [16] 



^ = cl^f + cU^f + cUg^f + 2cl^\^ + 2c\^af + 2^qf\f 

° l = 4(?i (3) ) 2 + 4 5 (^ (5) ) 2 + 2^ 5g ; (3) g; (5) + kWPf (20) 



with c" K denoting the collision matrix elements. Coefficients a r , a r , indicate the di- 
mensionless component of the electronic conductivity, the thermoelectric coefficient and the 



electric (a = e) or ion (a = i) thermal conductivity, respectively. From Eqs (16), (17) and 



(19), we get 

a a = <73(j£,x,t) + of(j£,x,t) ; a = (e,i) (21) 

Hence, also for collisional tokamak-plasma, the barycentric velocity does not appear explic- 
itly in the expression of the entropy production (see the footnote [42]). If the plasma becomes 
unstable and turbulent, this description becomes incomplete and supplementary contribu- 
tions to the expression of the entropy production should be added [26], [27]. However, we 
draw the attention to the fact that, also in the turbulent theory proposed in Ref. [10], the 
term which is supposed to be responsible of the intrinsic rotation is not linked to the velocity 
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flow but it depends only on its derivatives. In the work of Y. Kosuga et al, the term of 
the entropy source strength, a T , which is supposed to be responsible of the intrinsic plasma 
rotation can be brought into the bilinear ("flux x force") form : 

~t n 9 -< uii y 

a T = — =- -< u r uii > — !! — 22 

Here, r is the radial coordinate of the tokamak. In the mean field approximation we have 

^u r u|| y= - x / ^ - + ny (23) 

Here, x<t> represent the flow diffusivity and n^ d - the residual stress piece, respectively 
This latter contribution is not directly proportional to the toroidal rotation u^ or its radial 
gradient, but it is driven by the gradients of temperature, pressure and number density. 
Hence 



n 



7& 

° = ~2 

V thi L 



dr J r|1 dr 
According to the Diamond et al. theory, the contribution 

,d -< uii y 



Xd 9 ^) 2 -^ 9 ^' (24) 



= '^fi " ^ 



is a negative definite quantity and it is supposed to be responsible of the onset of the intrinsic 
"spin-up" of the plasma [27]. 

Our analysis leads then to the following conclusion: The global motion of tokamak-plasmas 
with barycentric velocity constitutes an intrinsic reversible processes. The thermodynamic 



theorems discussed in the Subsections II A and II B also apply to plasmas, magnetically 
confined in toroidal geometry. 

• Onset of the "Spin-up" of the Plasma from Rest 



We have shown that, if the Prigogine's hypothesis (see the Subsection II A) is fulfilled (and 
certain conditions satisfied), the results reported in Section [H] establish that a system may 
either stay at rest or undergo a global barycentric rigid motion without dissipation. Indeed, 
both situations are equally acceptable from the dynamical and thermodynamical point of 
view, because they satisfy Eq. (jl| and they perform a reversible process (i.e., they are not 
accompanied by dissipation). It will be the initial condition to select one between these two 
possibilities. However, many effects can drive a stationary plasma to rotate. The system may 
be subject, for example, to turbulence that can produce local Raynolds stress. The residual 
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stress, driven by a certain radial inhomogeneity, leads to the formation of the intrinsic 
rotation working in concert with the boundary conditions j2H]- Or, the excitation of modes 
driven by the plasma pressure gradient, may generate the inflow of angular momentum, 
whose source is at the edge of the plasma column [12]. Ultimately, internal hydrodynamic 
fluctuations, with high value of the wave number k, may also be enhanced by the nonlinear 
terms present in the balance equations and by the toroidal geometry [20j, [29]. Although, all 
of these mechanisms may explain the onset of the spin- up of the plasma from rest, we still 
have to understand the reason for which the plasma does not damp its rotation subsequently. 
Our task is, then, to estimate the solution of Eq. ([I]) with non vanishing initial flow velocity. 
Due to the predominance of the ion mass with respect to the electron mass, the velocity and 
momentum of the single fluid can be considered equal to the velocity and momentum of the 
ion fluid. Hence, we should solve the following evolution equation for the ion velocity 
duj 



171471;—— 

dt 



-Vpi - V ■ Hi + nZq e (e + - Ui x b) + R e + f ext (26) 



with initial condition Uoi 7^ 0. Here, c, B and E denote the speed of light, the magnetic, and 
the electric fields, respectively. R e and i ex t are respectively the interaction with the electron 
fluid and the external forces, per unit volume (other than electromagnetic). Since our aim 
is to study the intrinsic plasma rotation, we set i ext = 0. The scalar pressure p is defined 
through the equation of state Pi = n{Fi ~ nT and 7Tj denotes the ion viscous stress tensor. 
• Toroidal Rotation 

Our aim is to determine the toroidal angular frequency by applying the result shown at 
the end of the Subsection III Al and the Subsection III Bl In absence of external forces and 
by imposing the validity of the balance equation for the internal forces, the perpendicular 



velocity uj_ is immediately obtained from Eq. (26) by taking the vectorial product with B 



(see also the footnote [10]). At the leading order of the drift parameter, the perpendicular 
component of the flow equals the E x B drift and the diamagnetic flow (see, for example, 
Refs [3D] and [HI]) 

E x b b x VP . , , B 

U± " U ^^- + ^A Wlth h ^B (27) 

Here f2j denotes the ion gyrofrequency. Hence, at the lowest order of e and in the absence 
of poloidal flow, the toroidal rotation frequency, a^ n . c ., is 

d$ 1 dp 

^ n ' c ' dip nq e dip 

12 



where $ and ip denote the electrostatic potential and the poloidal flux, respectively. Hence, 
a driven toroidal plasma rotation is established in plasmas due to the formation of a radial 
electric field (in the presence of an externally applied magnetic field) and to a pressure 
gradient in the plasma. The plasma may then freely rotate in the toroidal direction, with 
frequency approximatively equal to 

Utf, = u^ n . c . + v^/R (29) 

with v denoting an arbitrary " departure" of the ion velocity from the barycentric velocity u, 
and R = Ro + r cos 6 (6 is the poloidal angle). In other words, plasma can freely rotate with 
an angular toroidal frequency, which may be increased with an arbitrary contribution Au^ = 



V(f,/R with respect to the neoclassical expression reported in Eq. (28). The initial conditions 
determine the value of v. In this respect, we mention that experimental measurements 
of toroidal rotation have been conducted in the L-mode, TCV ("tokamak a configuration 
variable") -plasmas, without external momentum injections [22]. It has been found that 
the neoclassical theory of toroidal rotation disagrees with the experimental observations 
qualitatively and quantitatively. It should be reminded, however, that if viscosity is taken 
into account the plasma may rotate toroidally with angular velocity approximatively given 



by Eq. (28) 



• Poloidal Rotation 

Let us now determine the poloidal angular frequency by applying the result discussed in the 



Subsection |II A| To determine the poloidal velocity flow, we should take into account the 
following two observations : 

j) In mechanical equilibrium, the barycentric velocity for plasma with no viscosity, inertial, 
nor external forces, obeys to the equation 

Vp ~ Vpi = nZq e (E + -3 x b) (30) 



By rewriting Eq. (30) in terms of the toroidal and the poloidal velocity components 



(u r ,u^,ue), we obtain the following relation between the angular frequencies 

B e R E r Vp-e r 

Ud = ^R~7 + TfT ~ nr7a Ft (31) 

tStf, r rJJ^ nr/jqeB^ 



with e r denoting the unit vector in the radial direction. Eq. (31) relates the poloidal and 
toroidal rotation to each other. In general it is not possible for the plasma to rotate poloidally 
without also rotating toroidally. 
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jj) In Eq. (26) we made use of the fact that all terms in the viscosity tensor are small 
(by a factor e) with respect to the other terms in the momentum equation, V ■ 7r Vp. 
However, this does not means that viscosity may always be neglected. In tokamaks, the 
parallel pressure gradient vanishes at the leading order and parallel viscosity plays a crucial 
role in determining the parallel barycentric velocity. 

An estimation of the order of magnitude of the barycentric velocity in the poloidal direction 
may be easily obtained by analyzing plasmas in the two-dimensional Velikhov model [19]. 
The model is constructed by taking into account considerations j) and jj). We consider a 
plane layer of plasma, which is in motion along the x-axis [which, according to j), simulates 
the presence of a "toroidal rotation"], but with no initial sliding motion in the y direction 
(this direction simulates a "poloidal rotation"). We suppose that the layer is subject to a 
magnetic field directed along the x-axis. In this case, the ion velocity, in the y-direction, sat- 
isfies the following simple equation, where the viscosity pressor tensor is no longer neglected 
[according to jj)] 



Eq. (32) reads 

(9 d\ d {1) du iy d {3) du ix 

The viscosity coefficients for ions are rj^ = 3njTj/(10fifrj) and rj^ = njTj/(20j) (with r» 
denoting the ion-ion collision time) [33] , [30] • In the low-collisional regime we have fijTj ^> 1 
and t)i can be neglected in this transport regime. Using this approximation, we obtain the 
following equation of continuity for the ion momentum in the y-direction 

dpiy d ( 1 d niTi\dui T\ d 2 Uj 



dt + Yx^ xVty = V iv ~n i Yx^m l )^x~~m l !w [ ' 

where Pi y = rriiUiy. We should now introduce the definition of the barycentric momentum in 
the y-direction, p y b. We have already mentioned that, for a system made by several species, 
the viscous tensor depends on the velocity of the reference frame [21]. In the single fluid 
description, the barycentric momentum is defined so as to satisfy the continuity equation 
where the transport coefficients of the viscous stress tensor are considered as constants 
and the compressibility of the fluid is negligible [T3]. Somehow, this definition "simulates" 
a global motion of the system. The viscous stress tensor is then constructed with the 
derivatives of the departure of the ion velocity from this barycentric velocity. Hence, in our 
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two-dimensional model, the dynamic equation for the barycentric momentum of the fluid in 
the y-direction reads 



'dp yb d 
,(3) 



-v 



(3) 



dx' 



(35) 



with f/ 3 ) = nT/{2VLi) ~ n;T;/(2^) = ^ o; . By imposing that the momentum of the single 
fluid is equal to the momentum of the ion fluid, we find that the system performs a "global" 
sliding motion ("rotation") for n ~ m and u y b = Ui y , with u y b given by 

0.5 &Ti 



u 



1 d riiTi 



yb 



(36) 



rii dx 2iriiVti rriiVti dx 

In toroidal coordinates (r, #,</>), a more accurate calculation for tokamak-plasmas with cir- 
cular cross section and large aspect ratio, gets (see, for example, [30] and |JT]) 

1.17 dTf 



U0 n .( 



mMi dr 



(37) 



In order of magnitude, we obtain the following neoclassical expression for the poloidal an- 
gular frequency 

1-17 dTi 

uen.c. ^ TV^T ( 38 ) 

rrriiili dr 

Notice that, with respect to the original work [TTJ], we have removed the ad hoc hypothesis 



stating that T/VLi is an adiabatic invariant. According to the results in the Subsection II A 



the plasma freely rotates poloidally, with angular velocity given by Eq. (38 ), without producing 



dissipation. Expression (38) roughly coincides with the experimental observations. However, 
since u y is the drift velocity of ions, according to Rosenbluth, a plasma moving with such 
velocity must be only marginally stable [31]. Taking into account Eq. (29) and Eq. (31 ), the 



toroidal rotation will tend to increment the neoclassical expression for the poloidal rotation 



[i.e., Eq. (38)], with an additional contribution equal to 



Bg V^ 

b2T 



(39) 



Therefore, according to the equation for plasmas in mechanical equilibrium, the toroidal 
rotation may cause the plasma to rotate poloidally at a speed faster than the neoclassical 



expression given by Eq. (38). However, the result discussed in the Subsection II A tell us that 
the departure of the poloidal angular frequency from the neoclassical value affects the value of 
the entropy production by contributing to dissipation. In addition, collisions tend to restore 
the neoclassical value. Hence, the magnitude of the poloidal angular frequency is determined 
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by the toroidal rotation (which drives the poloidal rotation) and by the neoclassical poloidal 
flow damping. Nevertheless, a variation of the poloidal angular rotation modifies, in its turn, 



the value of the toroidal angular frequency [through Eq. (31)], and so on. At the end, this 
mutual interplay breaks off when the toroidal and poloidal angular frequencies attains the 
values that minimize dissipation. 

IV. CONCLUSIONS 

From the results reported in Refs [IS] . [17j . and the Prigogine theorems, suitably readapted 
to toroidally confined dissipative plasmas, we showed that the global barycentric rotation 
does not contribute to dissipation. In case of negligible viscosity, we demonstrate that the 
global motion of tokamak-plasmas, satisfying the balance equation for the internal forces, 
is determined up to an arbitrary flow velocity. The validity of these results do not require 
that the system is at the mechanical equilibrium. In these conditions, the plasma may freely 
rotate in the toroidal direction with an angular frequency that may be higher than the neo- 
classical estimation, with an additional contribution equal to Au^ = v^/R (with v^ denoting 
an arbitrary velocity). This effect may help to explain the origin of the observed discrepancy 
between the neoclassical predictions and the experimental measurements in L-mode TCV 
plasmas. Plasma may also undergo a self-generated rotation in the poloidal direction with 
the angular frequency found by the neoclassical theory. Also in this case, the phenomenon 
is not accompanied by dissipation. However, from the Rosenbluth theory, a plasma moving 
close to this drift velocity is only marginally stable. The toroidal rotation may, then, force 
the plasma to increment the poloidal angular frequency with an additional contribution 



equal to Aug = Au^RBg/ '(rB^). The poloidal rotation, in its turn, will affect, via Eq. (31), 
the value of the toroidal angular frequency and so on. The eventual configuration will be 
attained when the toroidal and poloidal angular frequencies reaches the values that minimize 
dissipation. 

The next task should be, then, to determine the evolution equations for toroidal and poloidal 
angular frequencies and to analyze the stability of the mechanical equilibrium (or of the 
steady state solution). The study of this problem is beyond the scope of the present work. We 
have seen that, according to our description, two configurations of the plasma are compatible 
with the dynamical laws and thermodynamics. The plasma may either stay at rest or to 
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perform a global motion at the barycentric velocity. Both configurations are allowed by 
the thermodynamical principles because they are both related to reversible processes. As 
mentioned in the introduction, one of the leading theory suggests that the residual stress may 
spin- up the plasma from rest, acting in synergy with radial boundary conditions [36] (see also 
the footnote 02] )■ However, even if the above mentioned theory may suggest a mechanism 
able to interpret the origin of the rotation of the the tokamak-plasma in absence of external 
torque, it remains an open question the reason for which the plasma is able to self-sustain 
a global motion with a well determined angular frequency. The latter question is addressed 
in this work. We tackled this intricate problem by means of Prigogine's thermodynamic 
theorem, and by identifying the collisions as the only source of irreversibility or dissipation. 
Prigogine's theorem suggests that, once established the onset of the spin-up of the plasma, 
the system can self-maintain a global motion at the barycentric velocity since this process 
is not accompanied by dissipation. This theorem is a kind of thermodynamic equivalency 
version of the (mechanical) Galilean inertial principle, according to the scheme illustrated 
in FIG. [TJ We have remarked that the turbulent contributions, estimated by the Diamond 

Thermodynamics 
. , „ . . , < > Reversible Process 

Inertial Principle . ^ ... . 

% p Galileian Invanance 

Mechanics 

FIG. 1: In thermodynamics, a reversible process is "equivalent" to the Galileian invariance valid 
in mechanics, in the sense that both express the concept of inertial principle. 

et Al. theory, do not affect the validity of the results reported in Section [TTJ Moreover, 
the thermodynamic interpretation is perfectly in line with the experiment conducted on the 
DIII-D plasma, subject to external NBI [T8] . 

We conclude this section by mentioning that the paper cited in the footnote [43J, approaches 
the problem of obtaining a number of bounds of rotational levels in non spherical nuclei in 
a spirit similar to the one outlined above. 
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